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1.  BACKGROUND 


In  previous  work,  Segletes  (1991a.  1991b)  examined  entropy  changes  along  various  thermodynamic 
paths  of  simple  single-phased  materials  and  established  criteria  governing  three  modes  of  thermodynamic 
stability  of  the  Mie-Griineisen  equation  of  state  (EOS).  Since  these  criteria  are  used  as  a  starting  point 
for  the  current  analysis,  they  will  be  briefly  summarized  here. 

The  Mode  I  criterion  codifies  the  basic  notion  that  any  point  on  the  Hugoniot  in  p,v  space  must  rise 
more  steeply  (with  decreasing  v)  than  the  Rayleigh  line  connecting  the  Hugoniot  "foot"  to  that  point. 
Expressed  in  the  more  familiar  p.  p  space,  where  p  is  the  compression,  and  p  =  p/p^  -  1,  the  Mode  I 
criterion  becomes: 

dph/dp  ^  I 

Ph  pTTTp) 

where  P),  is  the  Hugoniot  pressure,  expressed  as  a  function  of  compression,  and  the  Hugoniot  "foot"  is 
taken  to  be  p^  =  0. 

The  Mode  II  criterion  makes  use  of  the  fact  that  thermodynamic  stability  requires  that  the  isentrope 
through  an  arbitrary  point  on  a  Hugoruot  curve  lie  between  the  Hugoniot  curve  itself  and  the  Rayleigh  line 
originating  from  the  Hugoniot  "foot"  (see  Figure  1).  As  a  result  of  enforcing  this  condition  upon  the 
Mie-Griineisen  EOS.  a  stability  relation  was  established  which  restricts  the  permissible  values  of  the 
GriiiKisen  coefficient  to  fall  in  the  range 

0  <  r  <  1  . 

P 

The  Mode  III  criterion  provides  an  inequality  involving  Griineisen  terms,  applicable  off  as  well  as  on 
the  Hugoniot,  which,  if  satisfied,  guarantees  that  the  slope  of  the  isentrope  through  that  point  is  of  the 
right  sign  (i.e.,  pressure  increases  with  increasing  compression).  It  is  given  by: 


r-t-1  ^  _i.iE 
1  +  p  r  dp 


(p-Ph)  + 


dPh 

dp 


(r/2) 

1  +p 


(Ph-Po)>o  • 


1 


Figure  1.  The  Mode  II  stability  criterion  is  based  upon  the  realization  that  the  slope  of  an  isentrooe 
through  a  point  must  lie  between  that  of  a  Hugoniot  through  the  tx)int  and  that  of  the  Ravleigh 


line  through  the  point  and  the  Hugoniot  "foot" 
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These  criteria  establish  new  guidelines  for  EOS  implementation  in  numerical  "hydrocodes. "  These 
hydrocodes  often  use  the  Mie-Griineisen  EOS  to  model  high  pressure  impacts,  involving  large 
deformations  and,  as  such,  can  be  sensitive  to  the  choice  of  functional  form  on  the  Griineisen  relation. 
Segletes  (1991a,  1991b)  shows  that  polynomial  Hugoniot  fits  in  the  literature  (Kohn  1969)  are  prone  to 
Mode  I  stability  violations  and  that  the  Griineisen  relations  existing  in  many  of  the  popular  hydrocodes 
have  generally  been  of  an  ad  hoc  nature  and  are  susceptible  to  Modes  11  and  111  violation  (e.g.,  EPIC 
[Johnson  et  al.  1978;  Johnson  &  Stryk  19921,  DYNA  [Whirley  et  al.  1992],  JOY  [Couch  1983),  MESA 
[Cagliostro  et  al.  1990],  and  to  a  lesser  extent,  HULL  [Matuska  &  Osborn  1987]). 

An  alternative  Griineisen  relation  was  proposed  by  Segletes; 


The  term  p  is  a  constant  parameter  which  must  be  chosen  greater  than  or  equal  to 

ro-2/p. 


where  r^,  is  the  Griineisen  value  at  the  Hugoniot  origin,  and  is  the  terminal  value  of  compression  (that 
compression  where  the  Hugoniot  pressure  becomes  unbounded).  This  alternative  form  is  guaranteed  to 
satisfy  the  Mode  II  criterion  at  all  valid  compression  states  and  was  found  to  greatly  diminish 
susceptibility  to  Mode  III  violation. 


For  general  materials,  these  rules  do  not  necessarily  hold  in  all  locations  in  p-v  thermodynamic  space. 
First  of  all,  because  the  current  discussion  is  limited  to  the  classical  case  where  the  Griineisen  coefficient 
is  taken  as  a  functional  of  volume  only.  Furthermore,  phase  changes  and  the  like  can  cause  unusual  local 
behavior,  for  example,  tensile  shocks  and  Mode  I  violations.  However,  these  unusual  regions  constitute 
the  exception  and  not  the  rule  over  the  full  thermodynamic  range  of  material  behavior.  In  a  sense,  then, 
the  rules  developed  are  applicable,  for  Griineisen  materials,  in  a  piece-wise  linear  fashion,  valid  between 
anomalous  thermodynamic  state  regions.  Furthetmore.  numerical  implementations  of  equations  of  state 
in  wave  propagation  codes  are  almost  universally  Griineisen  type  and  typically  continuous  and 
continuously  differentiable,  so  that  anomalies  of  the  type  discussed  (e.g.,  phase  changes)  are  not  modeled 
anyway.  Thus,  these  stability  criteria  take  on  added  importance  when  considering  numerical  aspects  of 
"production"  equation  of  state  modeling,  where  characterization  of  material  behavior  is  required  over  a 
wide  range  of  pressures  and  volumes. 


2.  THEORY 


In  the  current  effort,  a  fundamental  thermodynamic  rule  is  again  examined  to  see  what  influence,  if 
any,  is  exerted  upon  the  Mie-Griineisen  EOS.  The  rule  examined  is  simply  that  the  slope  of  a  Rayleigh 
line  -(dp/dv)  1 between  the  Hugoniot  origin  and  any  point  on  the  Hugoniot,  must  be  positive,  with 
increasing  compression.  Put  aitother  way,  the  rule  implies  that  the  pressure  at  a  shocked  state  must  exceed 
the  preshocked  pressure.  This  analysis  generalizes  the  Mode  III  criterion,  in  that  all  valid  preshock  states 
will  be  considered,  rather  than  just  those  arising  from  the  inflnitesimal  shock  (i.e.,  isentnopic)  condition. 

In  typical  numerical  implementations  of  the  Mie-Griineisen  EOS,  arbitrary  thermodynamic  states  of 
a  material  are  related  back  to  those  states  along  a  specified  reference  curve  at  the  same  specific  volume; 

P-Pref--^  (E-Eref)  • 

where  p  and  E  are  pressure  and  specific  internal  energy,  p„f  and  E^^f  are  the  pressure  and  specific  energy 
along  the  reference  curve,  and  y  =  v/T,  where  v  is  the  specific  volume.  In  hydrocodes,  the  reference 
curve  is  typically  chosen  as  the  Hugoniot  curve,  so  that  and  Ef^f  become  P|,  and  E^,  respectively.  In 
addition  to  the  Mie-Griineisen  relation,  the  following  discussions  will  also  make  use  of  the  shock  energy 
equation,  which  relates  the  specific  internal  energies  before  and  after  a  shock  front  to  the  pressures  and 
specific  volumes: 

^2  “  Ej  *  ^(Pi  +P2)(Vi  -Vj)  . 

Consider  a  shock  from  an  arbitrary  themtodynamic  state  "1"  to  state  "2,"  as  described  in  Figure  2. 
Note  that  the  states  corresponding  to  (with  the  same  specific  volume  as)  "1”  and  "2,"  but  along  the 
reference  Hugoniot  curve,  are  given  as  states  "hj"  and  "hj."  Also,  the  "foot"  of  the  reference  Hugoniot, 
which  is  almost  universally  the  ambient  state,  is  given  as  state  ”o."  Of  course,  Vj  >  V2  must  hold  to 
permit  shock  waves  into  simple  solids,  and  v^  >  Vj  must  hold  in  order  to  appropriately  apply  the 
Mie-Griineisen  EOS  with  Hugoniot  reference.  Thus,  the  relation,  v^  >  Vj  >  Vj,  will  hold  throughout  this 
analysis. 

For  the  current  analysis,  the  reference  state  ”o”  will  be  assumed  to  be  the  ambient  condition,  so  that 
pressure  and  internal  energy,  p^,  and  E^,,  will  both  be  set  to  zero.  The  ambient  state  is  almost  universally 
chosen  as  the  "foot"  of  experimental  Hugoniots,  so  that  this  choice  of  the  Hugoniot  reference  state 
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V 


Figure  2.  A  depiction  of  the  material  states  of  i  TCst:  state  "o"  is  the  reference  Hugoniot  "foot":  states 
"1"  and  "2"  are  the  ore-  and  post-shocked  states  under  consideration:  states  "hj"  and  "h^"  are 
those  states  on  the  reference  Hugoniot  corresponding  to  the  specific  volumes  of  states  "1"  and 
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provides  the  most  fruitful  assumption  for  study.  The  following  five  equations  will  now  be  considered; 
the  Mie-Griineisen  EOS  at  states  "1"  and  "2,"  and  the  shock  energy  equation  between  states  "o "  and  ”1, " 
"o"  and  "2,''  and  between  "1"  and  ”2,"  respectively; 


Pi  'Ph, -W 

Eh,*yPh,(''o-Vi) 

Ej-Ei  =I(Pi  +P2)(Vi  -Vj). 


It  is  desired  to  combine  these  equations  in  a  way  which  eliminates  all  energy  terms  (Ep  E2.  E(,^.  and  E,,^). 
Such  a  combination  may  be  obtained; 


P2 


V2-3  ('^i  -''2) 


Ph, 


V2  - 


*  Pi 


■''2^ 


-Ph, 


Vl  -  y(Vo  -  Vl) 


If  this  equation  is  used  to  express  pj  in  terms  of  p2,  then  the  positive  Rayleigh  slope  condition  being 
considered  (namely,  (P2  -  pj)  >  0),  is  given  as 


Ph, 

V2-4(''o-''2) 

-Ph, 

1 

0 

I 

P2 

[¥i  -¥2+(Vi  -V2)] 

¥1  +y(Vl  -V2) 
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The  denominator  oi  this  expression  is  always  positive  and  may  thus  be  eliminated  from  the  inequality 
without  changing  its  sense.  This  resulting  expression  may  be  used  to  isolate  the  terms  (p^  -  P2 )  in  the 
following  form: 


r  1  1 

(Phj-P2)(Vi-¥2+Vi-V2)  <  (Phj-Ph,)  ¥1  -  -jCVo-Vi)  +  Phj 


2  (''.-'’2) 


Though  this  paper  will  continue  to  derive  its  results  in  the  convenient  yff,  v  space,  results  may,  as  an 
alternative,  be  expressed  in  the  mote  familiar  F,  p  space  as  follows: 


(Phj  -  Pi)  ti'2  0  ^  Pi)  -  r,  (1  +  Pi)  -  r,  r2(P2  -  pi)]  <  (ph^  -  Ph,)  ^2  0  +  P2)  ) 

Ph, 

+  Fj  F2  (P2  -  Pi )  • 


r'lPi 
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This  inequality  (in  either  form)  defines  the  p.  v  (or  p,  p)  thermodynamic  space  where  p2  >  Pi 
true.  Those  state  pairs  "1"  and  "2"  which  fail  this  metric  imply  that  the  thermodynamic  path  of  shocking 
from  state  "1"  to  "2"  is  not  thermodynamically  admissible.  Considering  the  right-hand  side  of  the 
inequality,  it  is  clear  that  terms  like  pj,^,  (pj,^  -  Pj,^)  and  (Vj  -  Wj)  are  all  positive,  for  a  thermodynamically 
stable  condition.  The  remaining  term, 

¥i-j(Vo-Vi)  . 


when  converted  to  F,  p  space,  is  expressed  as 


1 

2 


From  the  Mode  11  criterion,  the  term  must  be  positive.  Since  F  and  v  are  always  positive,  this  term  and, 
in  fact,  the  complete  right-hand  side  of  the  inequality  will  always  be  positive. 


The  first  term  on  the  left-hand  side  of  the  stability  metric,  (pj^  -  pj).  will  be  positive  for  pressure 
states  below  the  reference  Hugoniot  and  negative  for  pressure  states  above  the  reference  Hugoniot.  The 
second  term  on  the  left,  (Vj  -  yj  +  a  function  of  specific  volume  only.  If  it  is  negative,  then, 

for  any  fixed  specific  volume,  pressure  states  far  enough  above  the  reference  Hugoniot  will  eventually 
violate  this  stability  metric.  If,  on  the  other  hand,  this  term  is  positive  for  a  given  specific  volume,  then 
pressure  states  far  enough  below  the  reference  Hugoniot  will  eventually  violate  the  metric.  The 
consequences  seem  startling,  since  they  imply,  unless  the  second  left-hand  term  is  precisely  zero,  that  the 
Mie-GrUneisen  EOS  will  necessarily  violate  Rayleigh  slope  stability  considerations  in  some  pressure  range. 
In  reality,  however,  valid  thermodynamic  states  are  bounded  from  below  by  the  "cold  curve";  that  p,  v 
curve  associated  with  the  absolute  zero  isotherm,  below  which  physically  plausible  states  do  not  exist. 
There  are  no  such  physical  upper  limitations  on  p,  v  states.  Therefore,  it  becomes  clear  that  the  only 
possible  state  of  the  stability  metric  which  can  make  physical  sense  is  to  have  the  second  left-hand  term 
be  nonnegative,  so  that  any  inadmissable  paths  lie  below  the  reference  curve — indeed,  far  enough  below 
the  reference  Hugoniot,  it  is  hoped,  so  as  to  lie  below  the  cold  curve. 

3.  MODE  IV  STABILITY  CRITERION 

This  stability  metric  forms  the  basis  for  a  new  constraint  on  the  Griineisen  relationship.  Adopting  the 
nomenclature  of  Segletes  (1991a,  1991b),  this  relationship  might  be  added  to  his  existing  list  and  thus  be 
termed  the  Mode  IV  criterion  for  EOS  stability: 


(¥i  -  ¥2)  +  (v,  -  Vj)  >  0  . 


Expressed  in  T,  p  space,  the  Mode  IV  criterion  becomes: 

^2(1  +  112)  -  r't(l  +  Mt)  +  Ti  r2  (P2  -  Ml)  ^  0  . 


Several  special  cases  of  this  criterion  are  examined  for  their  specific  functional  behavior.  For  example, 
if  the  prc-shocked  volume  is  the  ambient  specific  volume,  such  that  Vj  =  (or  correspondingly,  pj  =  0), 
then  the  resulting  constraint  on  the  Griineisen  relationship  becomes 


r  > 


1  +(i  +rj,)p 
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Note  that  this  criterion  provides  a  restriction  mote  constraining  than  that  of  the  Mode  II  criterion, 
0  <  r  <  2/p.  it  is  also  interesting  to  find  that  the  form  of  this  criterion  exactly  matches  the  functional 
form  of  the  Griineisen  relation  proposed  by  Scgletes  (1991a,  1991b).  namely: 


r 


To 

TTPiS’ 


where  p  is  a  constant  which  must,  in  the  present  case,  be  smaller  than  (1  +  in  order  to  satisfy  the 
inequality. 

As  a  result,  the  second  special  case  tested  does  not  fix  a  value  for  the  pre-shocked  specific  volume, 
v^,  but  rather  assumes  merely  that  the  functional  form  on  r  is  of  the  form,  involving  p,  proposed  above. 
The  result  in  this  case,  which  surprisingly  is  independent  of  both  prc-  and  post-shock  specific  volume, 
reveals  again  that 

P  ^  (I  +  r^) 

is  required  to  satisfy  the  Mode  IV  criterion.  On  the  other  hand,  the  Mode  II  criterion  shows,  when  this 
form  for  Griineisen  relation  is  adopted,  that  the  acceptable  range  of  P  is  given  by 

p.r.-K)  , 

2 

Recall  that  is  the  compression  at  which  the  Hugoniot  reference  function  becomes  unbounded.  Thus, 
to  adopt  the  form  of  Griineisen  relation  proposed  by  Segletes,  the  constraint  on  p,  by  combining  both  the 
Mode  IT  and  IV  criteria,  becomes 

Note  that  the  proposed  form  of  the  Griineisen  function,  with  P  so  constrained,  will  be  guaranteed  not  to 
violate  either  the  Mode  II  or  IV  criteria. 

The  final  special  case  investigated  for  the  Mode  IV  criterion  is  that  of  the  infinitesimal  shock,  where 
the  preshock  specific  volume,  Vj,  approaches  the  post-shock  condition,  \2-  special  case  turns  out  to 
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be  the  strictest  application  of  the  Mode  IV  criterion,  as  will  be  subsequently  shown,  and  is  thus  the  most 
useful  fonn.  The  differential  Mode  IV  criterion,  in  this  case,  becomes 


dv 

In  r.  p  space,  the  c*-'  ,ion  is  expressed  as 

dr  ^  r(i  +r) 

dp  (1  +p) 

As  dr/dp  identically  equals  the  right-hand  side  of  the  inequality,  Segletes’  form  can  be  recovered  through 
integration,  with  ^  equaling  (1  +  The  true  advantage  to  this  differential  form  of  the  Mode  IV 
criterion  is  that  any  proposed  Gtiineisen  function  may  be  evaluated  at  any  given  volume  v  (or  compression 
p).  rather  than  requiring  both  pre-  and  post-shocked  volumes  (v,  and  V2)  in  order  to  ascertain  stability, 
as  in  the  earlier,  nondifferential  form  of  the  criterion. 

It  was  stated  that  this  differential  form  of  the  criterion  is  the  strictest  form.  That  is  to  say,  it  fully 
encompasses  all  cases  covered  in  the  original  nondifferential  criterion.  To  prove  this  assertion,  one  must 
show  that  it  is  not  possible  for  a  functional  form  to  satisfy  the  dilTerential  criterion,  while  simultaneously 
failing  the  nondifferential  criterion.  It  may  seem  intuitive,  since  the  differential  form  constrains  the 
instantaneous  rate  of  change  of  the  Griineisen  function,  rather  than  merely  limiting  its  finite  change  over 
a  finite  specific  volume  interval.  However,  a  more  rigorous  proof  shall  be  employed  in  this  regard. 

One  may  devise  a  Griineisen  function  which  satisfies  the  differential  criterion  by  creating  a 
perturbation  from  that  criterion: 

^  +g(v);  g(v)  >  0  , 

dv 

where  g(v)  is  some  arbitrary,  nonnegative  perturbation  function.  Since  the  term  g(v)  is  always 
nonnegative,  this  functional  form  will  satisfy  the  differential  criterion  always.  The  solution  of  this 
differential  equation,  which  represents  (through  ip)  ail  plausible,  stable  Gruneisen  functions,  is: 

\|r  =  G(v)  +  Cj  -  V 

where  Cj  is  a  constant  of  integration  and  G(v)  is  the  integral  of  g(v)  with  respect  to  v.  Of  course,  this 
equation  is  subject  to  the  condition  that  (G(v)  +  Cj)  >  v,  in  order  to  keep  xp  (thus  F)  positive.  Since  g(v) 
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is  always  nonnegative,  G(v)  will  be  a  monotonically  increasing  function  of  v.  Substituting  this  vjt  function 
back  into  the  nondifferential  criterion  produces  the  result 

G(vi)  -  Gfvj)  S  0  . 

Since  G(v)  is  monotonically  increasing  and  >  Vj,  this  result  shows  that  the  nondifferendal  criterion  will 
always  be  satisfied  by  any  arbitrary  Griineisen  function  which  also  satisfies  the  differential  criterion. 
Therefore,  the  differential  criterion  is  also  a  rigorous  statement  of  Mode  IV  stability. 

Figure  3  depicts  the  consequences  of  the  Mode  IV  criterion,  by  illustrating  several  curves-of- 
maximum-descent.  These  curves  represent  the  steepest  possible  gradient  of  the  Griineisen  function,  short 
of  violating  the  Mode  IV  criterion. 

4.  ALTERNATIVE  MODE  III  STABILITY 

The  primary  consequence  of  Mode  IV  stability  is  that  any  and  all  unstable  thermodynamic  paths  are 
forced  to  occur  below  the  Hugoniot  reference  curve  of  the  Mie-Griineisen  EOS.  Figure  4  depicts  a  stable 
reference  Hugoniot  and  two  unstable  derivative  Hugoniots  that  can  occur  if  Mode  IV  stability  is  violated. 
The  stable  reference  Hugoniot  (for  aluminum)  is  based  upon  the  common  assumption  of  a  linear 
relationship  between  the  shock  velocity,  Uj,  and  the  particle  velocity  behind  the  shock,  Up.  The 
mathematical  form  is 


Us  =  +  s  Up  , 

and  the  Hugoniot  form  which  results  from  it  is  given  by 


PpCp^d  ^ 

[l-(S-l)p] 


Note  that  the  unstable  Hugoniots  lie  above  the  reference  Hugoniot,  in  thermodynamic  space  where 
Rayleigh  slope  violations  should  not  occur.  The  derivative  Hugoniots  were  generated  from  the  reference 
by  assuming  a  Griineisen  function  equal  to  r^/[l  +  (3  +  FjP  p],  which  was  shewn  in  the  previous  section 
to  be  clearly  unstable. 
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The  Mode  IV  criterion  defines  the  maximum  allov».tu.-'  ~"  /^-of-descent  of  the  GrOneisen 
function.  Several  such  curves  are  plotted  in  F.  u  space.  In  any  case,  the  Griineisen  value 
must  remain  below  a  value  of  (2/u).  as  dictated  by  the  Mode  11  criterion. 
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Figure  4.  A  depiction  of  unstable  derivative  Hugoniots  which  can  arise  from  a  Mode  IV  violation.  The 


stable  reference  Hugoniot  is  that  of  aluminum,  with 


a  value  for  r„  equal  to  2.09.  The  Griineisen  form  is  that  of  Segletes  (1991a.  1991b).  with  the 


value  for  B  set  to  (F,  +  3).  illustratively  chosen  to  lie  beyond  the  stable  limiting  value  of 


r.+  1 
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The  Mie-Griineiscn  EOS  form  under  consideration  is  defined  only  for  volumes  where  the  reference 
function  is  defined.  Since  the  Hugoniot  reference  function  is  only  defined  for  positive  compressions, 
discussions  throughout  this  section  are  limited  to  compressed  material  states  (at  specific  volumes  less  than 
the  ambient  specific  volume). 

Though  satisfying  the  Mode  IV  criterion  will  guarantee  that  thermodynamically  unstable  regions  will 
lie  below  the  reference  Hugoniot.  the  question  arises  as  to  how  far  below?  In  particular,  it  is  desired  to 
have  any  unstable  regions  lie  below  the  cold  curve  (i.e.,  the  zero-degree  isotherm),  where  they  would  pose 
no  problems.  Gases,  for  example,  have  their  zero-degree  isotherm  identically  at  p  =  0.  Solids,  on  the 
other  hand,  have  a  cold  curve  which  monotonically  increases  in  pressure  with  compression  and,  except 
for  a  small  region  near  ambient  density,  will  always  lie  above  p  =  0. 

Thus,  if  it  can  be  shown  that  all  thermodynamic  instabilities  lie  below  the  cold  curve,  while  satisfying 
the  Mode  IV  criterion,  then  one  may  conclude  that  such  instabilities  are  no  longer  relevant,  since  they  lie 
in  a  nonachievable  range  of  the  EOS  range. 

Recall  the  relation  required  to  keep  the  Rayleigh  slope  positive: 


Ph, 


V2- 


-  Ph, 


Vl-y(Vo'Vi) 


P2 


[vi -V2  +  V1 -Vj]  >  0  , 


where  V|  >  Vj.  The  differential  form  of  this  equation  follows  as: 


d 


v-yK-v) 


ll.p 

fii.il 

JJ 

[dv  J 

>  0  . 


The  limiting  case  for  the  stability  of  this  equation  is  when  p  is  at  a  minimum,  namely,  along  the  cold 
curve.  Pg.  Though  not  shown  here,  a  similar  analysis  may  be  employed  to  that  derived  previously,  which 
shows  the  differential  criterion  to  fully  encompass  the  nondifferential  criterion.  If  one  now  takes  y  times 
the  Mie  Griineisen  equation  of  state,  with  the  ambient  Hugoniot  as  its  reference,  where 
E^  =  1/2  p,,  (Vq  -  v) ,  one  obtains: 


14 


PV-E=Ph 


The  latest  two  terms  may  be  combined  by  talcing  the  derivative  of  this  latest  expression  and  substituting 
into  the  previous  expression,  eliminating  the  term 


d/dv 


The  result  is 


/ 


dvir  ^Pe 
dv  dv 


dE, 

dv 


dv 


>0 


Utilising  the  facts  that  P(,  =  -dE^./dv,  and  that  vjr  is  positive,  this  equation  ultimately  simplifies  to 


One  may  thus  conclude  that  a  monotonic  cold  pressure  curve  is  sufficient  to  keep  the  Rayleigh  slope 
stable. 

Recall  however,  that  the  differential  Rayleigh  slope  (i.e.,  for  infinitesimal  shocks)  coincides  with  the 
isentiopic  and  Hugoniot  slopes.  This  Rayleigh  stability  condition  is  thus  equivalent  to  an  isenttope 
stability  condition,  which  was  the  basis  of  Segletes’  Mode  III  criterion.  The  difference  lies  in  the  fact 
that  the  original  criterion  was  expressed,  not  in  terms  of  cold  curve  behavior,  but  rather  in  terms  of 
Hugoniot  behavior.  Nonetheless,  the  two  criteria  are  alternatives,  given  that  the  Hugoniot  and  cold  curve 
functions  may  be  derived  from  each  other  when  the  Griineisen  function  is  known. 


5.  SUMMARY 

The  Griineisen  parameter  within  the  Mie  Griineisen  EOS  is  often  treated  as  an  independent  function 
of  specific  volume,  judging  from  forms  employed  in  the  hydrocode  literature.  This  report  and  its 
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predecessor  (Segletes  1991a,  1991b)  strive  to  demonstrate  the  coupling  which  exists  among  EOS 
parameters  (including  the  Gruneisen  parameter),  and  in  so  doing,  set  forth  criteria  which  govern  various 
aspects  of  thermodynamic  stability.  The  types  of  EOS  behavior  which  may  be  expected  from  violations 
of  the  various  criteria  have  also  been  depicted. 

In  this  report,  an  inequality  was  established  which  quantifies  the  interrelation  of  thermodynamic 
quantities  necessary  to  prevent  the  occurrence  of  a  Rayleigh  slope  instability  (whereby  pre-shock  pressure 
erroneously  exceeds  post-shock  pressure).  This  inequality  was  expressed  in  the  numerically  convenient 
y,  V  space,  as  well  as  the  more  popular  T.  p  space  (see  Table  1). 

By  making  use  of  this  relation,  two  criteria  were  developed  (denoted  Mode  IV  and  alternate  Mode  III). 
Satisfaction  of  the  Mode  IV  criterion  will  guarantee  that  Rayleigh  slope  violations  can  only  possibly  occur 
below  (and  not  above)  the  reference  Hugoniot  curve  (see  Table  2).  Such  a  condition  is  necessary  if  there 
is  any  hope  of  keeping  stability  violations  below  the  thermodynamic  cold  curve  in  a  thermodynamically 
unachievable  range. 

The  alternate  Mode  III  critenon  establishes  a  relation  which  will  keep  differential  Rayleigh  slope 
instabilities  (thus,  isentropic  instabilities  below  the  cold  curve  (see  Table  3).  Simply  stated,  the  cold  curve 
must  monotonically  decrease  with  increasing  volume.  This  criterion  offers  an  alternate  form  over 
Seglotes’  original  Mode  III  criterion  which  expresses  isentropic  stability  in  terms  of  Hugoniot  rather  than 
cold  curve  behavior. 


Table  1.  Criterion  to  Avoid  Rayleigh  Slope  Instabilities 
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Table  2.  Mode  IV  Criterion  to  Keep  Possible  Rayleigh  Slope  Instabilities  Below  Hugoniot 


nondifferential  fonn  (y,  v  space) 

(\pj  -  V2)  +  (vi  -  V2)  >  0 

nondifferential  fonm  (F,  p  space) 

Fjfl  +  M2) 

-  Fjd  +  pj)  +  Fj  Fj  (P2  -  Pi)  ^  0 

differential  form  (vp,  v  space) 

ill  >  -1 

dv 

differential  form  (F,  p  space) 

dF  ^  _  F(1  +F) 
dp  (1+p) 

Table  3.  Alternate  Mode  III  Criterion  to  Eliminate  Rayleigh  and  Isentropic  Slope  Instabilities 
Above  Cold  Curve 


The  stability  criteria  resulting  from  this  work  provide  guidance  to  those  researchers  who  wish  to  make 
use  of  the  Mie-Gruncisen  EOS,  in  hydrocodes  or  other  media.  Adherence  to  the  Modes  III  and  IV  criteria 
remove  the  threat  of  Rayleigh  slope  instability  over  the  plausible  equation  of  state  domain.  Finally,  a 
better  understanding  of  the  interdependence  of  the  Mie-Griincisen  equation  of  state  parameters  may  be 
obtained  through  the  understanding  of  the  thermodynamic  stability  criteria  described  in  this  report. 


Intentionally  left  blank. 
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ATTN:  M.  T.  Gonzalez 

D.  Markel 

P.O.  Box  1590/14910  Ariline  Rd. 
Rosharon,TX  77583-1590 

2  Aerojet  Ordnance/Exp.  Tech.  Ctr. 
ATTN:  P,  Wolf 
G.  Padgett 
1100  Bulloch  Blvd. 

Socorro,  NM  87801 

1  Aerojet  Electro  Systems  Co. 

ATTN:  Warhead  Systems,  J.  Carleone 
P.O.  Box  296 
Azusa,  CA  91702 

2  Physics  International 
ATTN:  R.  Funston 

L.  Garnett 
2700  Merced  St. 

P.O.  Box  5010 

San  Leandro,  CA  94577 


2  Lockheed  Missile  &  Space  Co.,  Inc. 
ATTN:  S.  Kusumi  Q-Sl-H.  Bldg.  157) 
J.  Philips  (0-54-50) 

P.O.  Box  3504 
Surmyvale,  CA  94088 

1  Lockheed  Missile  &  Space  Co.,  Inc. 

ATTN:  R.  A.  Hoffman 
Santa  Cruz  Facility 
Empire  Grade  Rd. 

Santa  Cruz,  CA  95060 

1  Boeing  Corp. 

ATTN:  T.  M.  Murray  (MS-84-84) 

P.O.  Box  3999 
Seattle.  WA  98124 

2  Mason  &  Hanger  -  Silas  Mason  Co. 
ATTN:  T.  J.  Rowan 

C.  Vogt 

Iowa  Army  Ammunition  Plant 
Middletown.  lA  52638-9701 

I  Nuclear  Metals,  Inc. 

ATTN:  J.  Schreiber 
2229  Main  St 
Concord,  MA  01742 

1  Lockheed  Engineering  &  Space  Sciences 
ATTN:  E.  Cykowski  (MS  B-22) 

2400  NASA  Rd.  1 
Houston,  TX 

3  Dyna  East  Corp. 

ATTN:  P.  C.  Chou 

R.  Ciccarelli 
W.  Ris 
3201  Arch  St. 

Philadelphia,  PA  19104 

2  Southwest  Research  Inst. 

ATTN:  C.  Anderson 

A.  Wenzel 
6220  Culebra  Rd. 

P.O.  Drawer  28510 
San  Antonio,  TX  78284 
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Battelle  -  Columbus  Laboratories 
ATTN:  R.  Jameson 

S.  Golaski 
505  King  Ave. 

Columbus,  OH  43201 

AUiant  Techsystems.  Inc. 

ATTN:  G.  R.  Johnson,  MNn-2925 

T.  Holmquist,  MNll-2720 
R.  Stryk 

600  Second  Sl,  NE 
Hopkins,  MN  55343 

S-Cubed 

ATTN;  R.  Sedgwick 
P.O.  Box  1620 
La  Jolla,  CA  92038-1620 

California  Research  and  Technology,  Inc. 
ATTN:  R.  Fanzen 
D.  Orphal 
R.  Brown 
5117  Johnson  Dr. 

Pleasanton,  CA  94566 

California  Research  and  Technology,  Inc. 
ATTN:  M.  Majerus 
P.O.  Box  2229 
Princeton,  NJ  08543 

Orlando  Technology,  Inc. 

ATTN:  D.  Matuska 
J.  Osborn 
P.O.  Box  855 
ShaUmar,FL  32579 

Kaman  Sciences  Cotp. 

ATTN:  D.  Bameue 
D.  Elder 
P.  Bussel 
P.O.  Box  7463 
Colorado  Spring,  CO  80933 

LTV  Aerospace  and  Defense  Co. 

LTV  Missiles  and  Electronics  Group 
ATTN:  G.  L.  Jackson  (WT-71) 

P.O.  Box  650003 
Dallas,  TX  75265-0003 


No.  of 

Copies  Organization 

2  General  Research  Corp. 

ATTN:  A.  Charters 
T.  Menna 
P.O.  Box  6770 

Santa  Barbara,  CA  93160-6770 

1  General  Dynamics 

ATTN;  J.  H.  Cuadros 
P.O.  Box  50-800 
Mail  Zone  601-87 
Ontario,  CA  91761-1085 

1  Defense  Technology  International,  Inc. 
ATTN:  D.  E.  Ayer 

The  Stark  House 
22  Concord  SL 
Nashua,  NH  03060 

2  E.  I.  DuPont  De  Nemours  &  Co. 
ATTN:  B.  Scott 

L.  Minor 

Security  Director,  Legal  Dept 
P.O.  Box  1635 
WDmington,  DE  19899 

1  SRI  International 

ATTN:  Dr.  L.  Seaman 
333  Ravenswood  Ave. 

Menlo  Park,  CA  94025 

1  Martin  Marietta  Aerospace 

ATTN:  D.  R.  Bragg 
P.O.  Box  5837,  MP  109 
Orlando,  FL  32855 

1  Northrop  Corp. 

Electro-Mechanical  Div. 

ATTN:  D.  L.  Hall 
500  E.  Orangethorpe  Ave. 

Anaheim,  CA  92801 

3  Boeing  Aerospace  Co. 

Shock  Physics  &  Applied  Math 
Engineering  Technology 
ATTN;  R.  Helzer 
J.  Shrader 
T.  Murray 
P.O.  Box  3999 
SeatUe,  WA  98124 
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1  D.  R.  Kennedy  and  Associates  Inc. 

ATTN:  D.  Kennedy 
P.O.  Box  4003 
Mountain  View.  CA  94040 

1  Livermore  Software  Technology  Corp. 

ATTN:  J.  O.  Hallquist 

2876  Waverly  Way 
Livermore,  CA  94550 

2  University  of  Minnesota 

Army  High  Performance  Computing  Research 
Center 

ATTN:  G.  Sell 

D.  Austin 

1100  Washington  Ave.,  S. 

Minneapolis,  MN  55415 

3  U.S.  Army  Research  Laboratory 
ATTN:  T,  Chou 

C.  S.  White 
J.  McLaughlin 

Watertown,  MA  02172-0001 
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Aberdeen  Proving  Ground 

37  Dir.  USARL 

ATTN:  AMSRL-WT-NC, 

R.  Loterro 
K.  Opalka 

R.  Pearson 

S.  Schraml 
AMSRL-WT-TA, 

W.  Bruchey 
G.  Bulmash 
J.  Dehn 

G.  Filbey 

H.  Meyer 
AMSRL-WT-TB, 

R.  Frey 
W.  Lawrence 

J.  Starkenberg 
AMSRL-WT-TC, 

T.  Bjerke 
R.  Coates 

F.  Grace 

K.  Kimsey 

M.  Lampson 

D.  Scheffler 
B.  Sorensen 

R.  Summers 

E.  Walker 
W.  Walters 

ATTN:  AMSRL-WT-TD, 

A.  Dietrich 
K.  Frank 
J.  Harrison 

N.  Hufnngton 

G.  Randers-Pehrson 
J.  Santiago 

M.  Scheidler 

S.  Segletes  (3  copies) 

J.  Walter 

T.  Wright 
AMSRL-WT-W, 

A.  Prakash 
A.  Zielinski 
AMSRL-WT-P, 

G.  Gazonas 

Aberdeen  Proving  Ground,  MD  21005 
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Defense  Research  Establishment  Suffield 
ATTN:  C.  Weickert 
D.  MacKay 

Ralston,  Alberta,  TOJ  2NO  Ralston 
CANADA 

Defense  Research  Establishment  Valcartier 

ATTN:  N.  Gass 

P.O.  Box  8800 

Couicelette,  PQ,  GOA  IRO 

CANADA 

Canadian  Arsenals,  LTD 
ATTN:  P.  Pelletier 
S  Montee  des  Arsenaux 
Villie  de  Gardeur,  PQ,  J5Z2 
CANADA 

Ernst  Mach  Institute 
ATTN:  A.  J.  Stilp 
Eckerstrasse  4 
D-7800  Freiburg  i.  Br. 

GERMANY 

lABG 

ATTN:  H.  J.  Raatschen 
W.  Schittke 
F.  Scharppf 
Einsteinstrasse  20 
D-8012  Ottobrun  B.  MUnchen 
GERMANY 

Royal  Armament  R&D  Establishment 
ATTN:  I.  CulUs 
Ft.  Halstead 

Sevenoaks,  Kent  TN14  7BJ 
ENGLAND 

Centre  d’Etudes  de  Gramat 
ATTN:  SOLVE  Gward 
46S00  Gramat 
FRANCE 

Centre  d’Etudes  de  Vaujours 
ATTN:  PLOTARD  Jcan-Paul 
Boice  Postale  No.  7 
77181  Country 
FRANCE 


Copies  Organization 

1  PRB  S.A. 

ATTN:  M.  Vansnick 
Avenue  de  Tervueren  168,  Bte.  7 
Brussels,  B-llSO 
BELGIUM 

1  AB  Bofors/Ammunition  Div. 

ATTN:  J.  Hasslid 
BOX  900 
S-691  80  Bofors 
SWEDEN 

1  TNO  Prins  Maurits  Lab. 

ATTN:  H.  J.  Passman 
P.O.  Box  45 
2280  AA  Rijswijk 
Lange  Kleliweg  137 
THE  NETHERLANDS 

1  ISL 

ATTN:  CHANTERET  Pierre-Yves 
F68301  Saint-Louis 
C6dex,  12,  rue  de  I’lndusuie 
BP.  301 
FRANCE 

1  IFAM 

ATTN:  Lothar  Meyer 
Lesumo*  Heerstrasse  36 
Bremen  77 
GERMANY 


Intentionally  left  blank. 
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USER  EVALUATION  SHEET/CHANGE  OF  ADDRESS 

This  Laboratory  undertakes  a  continuing  effort  to  improve  the  quality  of  the  reports  it  publishes.  Your 
comments/answers  to  the  items/questions  below  will  aid  us  in  our  efforts. 

1.  ARL  Report  Number  ARL-TR-234 _ Date  of  Report  October  1993 _ 

2.  Date  Report  Received _ 

3.  Does  this  report  satisfy  a  need?  (Comment  on  purpose,  related  project,  or  other  area  of  interest  for 

which  the  report  will  be  used.) _ 


4.  Specifically,  how  is  the  report  being  used?  (Information  source,  design  data,  procedure,  source  of 
ideas,  etc.) _ 


5.  Has  the  information  in  this  report  led  to  any  quantitative  savings  as  far  as  man-hours  or  dollars  saved, 
operating  costs  avoided,  or  efficiencies  achieved,  etc?  If  so,  please  elaborate. _ 


6.  General  Comments.  What  do  you  think  should  be  changed  to  improve  future  reports?  (Indicate 
changes  to  organization,  technical  content,  format,  etc.) _ 


Organization 

CURRENT  Name 

ADDRESS  _ 

Street  or  P.O.  Box  No. 

City,  State,  Zip  Code 

7.  If  indicating  a  Change  of  Address  or  Address  Correction,  please  provide  the  Current  or  Correct  address 
above  and  the  Old  or  Incorrect  address  below. 

Organization 

OLD  Name 

ADDRESS  _ 

Street  or  P.O.  Box  No. 

City,  State,  Zip  Code 

(Remove  this  sheet,  fold  as  indicated,  tape  closed,  and  mail.) 

(DO  NOT  STAPLE) 


Department  of  the  army 


OFRCIAL  BUSINESS 


BUSINESS  REPLY  i\L\IL 

RUST  CUBS  BSMT  Ho  Ml,  >K.  MU 

Postage  i*ill  &e  paid  Oy  addiessee 


Director 

U.S.  Army  Research  Laboratory 
ATTN:  AMSRL-OP-CI-B  (Tech  Lib) 
Aberdeen  Proving  Ground.  MD  21005-5066 
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